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Summary. We present mid-inlrared interlerometric observations ol the Herbig Ae 
star R CrA obtained with MIDI at the VLTl using several projected baselines in 
the UT2-UT3 configuration. The observations show resolved circumstellar emission 
in the wavelength range 8-13 /^m on a ~6 AU scale with a non-symmetric intensity 
distribution, providing support for an inclined disk geometry. Visibilities are best 
fitted using a uniform ring model with outer radius in the range 6-10 AU, in the 
wavelength range 8-13 micron. The inclination of the ring with respect to the plane 
of the sky is found to be ~ 45° , consistent with the 40° suggested from near-infrared 
imaging polarimetry (Clark et al. 2000, MNRAS, 319, 337). 



1 Introduction 

The presence of circumstellar disks around the intermediate mass (M<5M0) 
Herbig Ae stars is supported by a large body of observational evidence [1][2]. 
More specifically, while the observed spectral energy distribution (SED) of 
such stars can be explained by both a disk-like distribution of material 
(e.g. [3]) and other geometries like envelopes (e.g. [4]), clear evidence for cir- 
cumstellar disks comes from resolved flattened structures observed by in- 
terferometry at millimeter, near-IR and recently also mid-IR wavelengths 
(e.g. [5][6][7]). 

R CrA is a bright (100 L©) young Herbig A5e star, located at the center 
of a small cluster (The Coronet [8] [9]) at 130 pc [10]. Several characteristics 
indicate the presence of a circumstellar disk around R CrA : a flat mid-IR 
to far-IR/mm SED'^ ([2]), a broad silicate emission feature [12], a UX Ori 
Type [13], a high degree (8%) of optical linear polarisation [14] , the possi- 
ble association with a extended molecular outflow [15], [16] as well as with 
several Herbig-Haro systems [17], and a near-infrared reflection nebulosity 
whose resolved spatial polarization is consistent with a bipolar outflow being 
truncated by an evacuated spherical cavity [20] . 



^although most of the mm-excess is actually from the nearby embedded infrared 
source IRS7 [11] and source confusion in the large IRAS beams might be an issue. 
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Fig. 1. Left panel: Wavelength averaged visibilities as a function of ^^^^radius and 
best-fit face-on uniform ring brightness distribution (Rout=6.3 AUIq'j). Dashed and 
dotted lines are respectively the inner and outer boundary values of Rout- Right 
panel: Contour plot of the best-fit inclined uniform ring model whose parameters 
are listed in Table 1. The contour levels are for visibilities of 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, 0.9. East is up, north is right. 

2 Observations and Calibration 

R CrA has been observed with MIDI at the VLTI between 2004 July 8 and 
30. We obtained 6 sets of spectrally dispersed (Prism, R=30) visibilities in 
the spectral range 8-13 /im with the 47 m long northeast UT2-UT3 base- 
line. Projected baseline lengths and PAs are 27.4 m, 46.6 m, 46.2 m, 45.6 m, 
40.1m, 39.6m and 63?6, 30?5, 9?4, 43?4, 55?6, 56?1, respectively. Each data 
set has been reduced with the MIA package [18]. The raw visibilities were 
calibrated by observing the source HD 173484, whose adopted angular size is 
3.35 ±0.31 mas (MIDI calibrator list). We evaluated the error on the visibili- 
ties by comparing the instrumental visibility obtained with different calibra- 
tors observed during the same night. This leads to a typical relative uncer- 
tainty of 15%, a value we adopt in the following. The acquisition images from 
both telescopes show that R CrA is unresolved at 8.7 /im, i.e. at a resolution 
limit of A/£'=0'.'27, corresponding to ~30AU at 130 pc. 

3 Evidence for an inclined disk structure 
3.1 First approach 

We compared the observed visibilities with those derived from a simple geo- 
metric disk model chosen to be a uniform ring brightness distribution. The 
inner radius was fixed to a dust sublimation radius of 0.3 AU, typical for a 
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Herbig AeBe star of that luminosity [19]. In a first step, only the wavelength 
averaged visibilities were fitted. We will see later that this rough approxima- 
tion can be justified by the relative flatness of the visibilities observed on R 
CrA in the 8-13/xm wavelength range. As a start, and for comparison pur- 
poses, we took a ring with zero inclination (face-on), so that the outer radius 
Rout was the only free parameter. 

Left panel of Fig. 1 shows the visiblities plotted as a function of M^^radius 
together with the best-fit model. For the latter, Rout=6.3 AU jlg ^ and the re- 
duced x^=2.4. Two conclusions can be drawn from this preliminary analysis. 
First, the mid-IR emission of R CrA is clearly resolved into a structure of 
6 AU radius. Second, this structure is asymmetric which hints at a inclined 
disk geometry. Including the inclination in the model leads to a better fit, with 
a reduced x^=l-l- The best-fit model is shown in contour plot in the right 
panel of Fig. 1 together with the ?7,?>coverage of our observations. The fitted 
parameters are Rotif=8.6 AU^^ g, an inclination of i=47° I22 with respect to 
the plane of the sky, and a semimajor axis position angle* PA=152° t^^. 

3.2 Disk structure and orientation 

Spectrally-dispersed visibilities as measured by MIDI are likely to provide 

additional information, because of the extension of ^^^^coverage. However, 
one should keep in mind that there is in principle a degeneracy between 
spatial brightness distribution (morphology) and dust spectral emission fea- 
tures. A fortunate observation, in this respect, is that the ISO-SWS spectra 
are rather featureless apart from a broad 9.7 /xm amorphous silicate feature, 
i.e. no prominent PAH feature at 8.6 /xm and no 11 /xm complex feature [12]. 
This allows us to derive meaningful results by fitting our simple uniform ring 
model to the observed spectrally-dispersed visibilities. Fig. 2 shows the set 
of observed spectrally-dispersed visibilities together with the best-fit face-on 
and inclined ring models. Given the possibility that the observed visibilities 
may be due to the presence of a close companion^, we additionally fit a bi- 
nary model to the complete data set. The fitted parameters as a function 
of wavelength are presented in Fig. 3, while their wavelength-averaged values 
can be found in Table 1. The quoted errors correspond to a variation of re- 
duced of unity. While we can rule out with a high degree of confidence a 
binary model, the best agreement with these mid-IR visibility data is found 
for an inclined uniform ring model. 

Fitting a uniform ring model to the total data set for each spectral channel 
as we did leads to a monotonic variation of Rout with wavelength (Fig. 3). 
Smaller outer radii are found at shorter wavelengths. This explains the fact 
that emission closer to the star arises from hotter dust and could in principle 
be used to derive the temperature profile of the circumstellar disk. 

"^oriented east of north. 

^R CrA is flagged in Hipparcos as "stochastic binary" and "suspected nonsingle" 
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Fig. 2. Spectrally-dispersed visibilities between with 8-13 /im (diamonds) with the 
best-fit models (binary, face-on and inclined uniform rings). 




Fig. 3. Best fit parameters as a function of wavelength for respectively a face-on 
(left panel) and inclined (right panel) uniform ring model. 
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Table 1. Best-fit models based on spectrally-dispersed visibilities. Uniform ring 

parameters arc wavelength-averaged values. The binary model separation, PA and 
brightness ratio are quoted in the Rout, PA and i columns, respectively. 







reduced 


Ro^t (AU) 


if) 


PA (°) 


Binaxy 

Face-on uniform 
Inclined uniform 


ring 
ring 


2.4 
2.2±0.7 
1.3±0.6 


q 1 +0.2 
-0.2 

8.48 ll;f. 


0.44 
44 If, 


Of! +10 

oo -10 

152114 



The inclination of the ring with respect to the plane of the sky is 44° 
This is fully eonsistent with the inclination of the plane (40°) perpendicular 
to the symmetry axis of the bipolar reflection nebula as suggested from NIR 
imaging polarimetry [20]. However, the derived position angle of the disk sym- 
metry axis (62° is only marginally consistent with R CrA as the driving 
source for the extended NE-SW molecular outflow at PA- 30° [15], [16]. The 
same conclusion holds for the direction of the closest HH objects with respect 
to R CrA [17] (HH96: ~38°, HH97W: ~34°, HH98: ~33°, HHIOO: ~34°, 
HH99A: ~51°, HH104A: ~95°). 

Acknowledgment: We are grateful to B. Wilking for his help. GM is supported 
by DFG grant ME 2061/3-1. 



References 

1. Natta, A., Grinin, V., Mannings, V., Protostars and Planets IV, p. 559 (2000) 

2. Hillenbrand, L.A., Strom, S.E., Vrba, F.J. et al., ApJ, 397, 613 (1992) 

3. DuUemond, CP., Dominik, C. and Natta, A., ApJ, 560, 957 (2001) 

4. Hartmann, L., Kenyon, S.J., Calvet, N., ApJ, 407, 219 (1993) 

5. Mannings, V. & Sargent, A.I., ApJ, 490, 792 (1997) 

6. Eisner, J. A., Lane, B.F., Hillenbrand, L.A., ApJ, 613, 1049 (2004) 

7. Leinert, Ch., Van Boekel, R., Waters, L.B.F.M. et al., A&A 423, 537 (2004) 

8. Taylor, K.N.R. & Storey, J.W.V., MNRAS, 209, 5 (1984) 

9. Wilking, B., McCaughrean, M., Burton, M.G. et al., AJ, 114, 5 (1997) 

10. de Zeeus, P.T., Hoogcrwerf, R., de Bruijne, J.H.J., AJ, 117, 354 (1999) 

11. Choi, M. & Tatematsu, K., ApJ, 600, L55 (2004) 

12. Acke, B. & Van den Ancker, M.E., A&A 426, 151 (2004) 

13. DuUemond, CP., Van den Ancker, M.E., Acke, B. et al., ApJ 594, L47 (2003) 

14. Bastien, P., ApJ, 317, 231(1987) 

15. Levreault, R.M., ApJSS, 67, 283 (1988) 

16. Anderson, I.M., Harju, J., Knee, L.B.G., A&A, 321, 575 (1997) 

17. Wang, H., Mundt, R., Henning, T. et al., ApJ, 617, 1191 (2004) 

18. Kohler, R. and Jaffe, W.: "MIA-hEWS, the software for MIDI data reduction". 

19. Monnier, J.D., Millan-Gabet, R., Billmeier, R. et al., ApJ, 624, 832 (2005) 

20. Clark, S., McCall, A., Chrysostomou, A. et al., MNRAS, 319, 337 (2000) 



